Quasiparticle properties of strongly correlated electron systems with itinerant 
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A brief account of the zero temperature magnetic response of a system of strongly correlated 
electrons in strong magnetic field is given in terms of its quasiparticle properties. The scenario is 
based on the paramagnetic phase of the half-filled Hubbard model, and the calculations are carried 
out with the dynamical mean field theory (DMFT) together with the numerical renormalization 
group (NRG). As well known, in a certain parameter regime one finds a magnetic susceptibility 
which increases with the field strength. Here, we analyze this metamagnetic response based on 
Fermi liquid parameters, which can be calculated within the DMFT-NRG procedure. The results 
indicate that the metamagnetic response can be driven by field-induced effective mass enhancement. 
However, also the contribution due to quasiparticle interactions can play a significant role. We put 
our results in context with experimental studies of itinerant metamagnetic materials. 

PACS numbers: 71.10.Fd, 71.27.+a,71..30.+h,75.20.-g, 71.10.Ay 



I. INTRODUCTION 

The interplay of strong correlation physics and mag- 
netic behavior in itinerant electronic systems has been a 
fascinating subject for many years. At low temperature it 
is often possible to describe the response of such systems 
in terms of the low energy excitations and quasiparticle 
properties such as in a Fermi liquid picture. The ratio of 
the spin susceptibility of the interacting system Xs and 
that of the non-interacting system x° is then given by 
the expression 



Xs _ m*/mo 



1 + F- 



(1) 



where m* /mo is the ratio of effective and bare electronic 
mass, and Fq is the lowest order asymmetric Landau pa- 
rameter, which accounts for quasiparticle interactions. A 
special kind of response is metamagnetism, which we de- 
fine here as the existence of a regime where the system's 
differential susceptibility Xs = dM/dH increases with 
magnetic field H, i.e. dxs/dH > 0, for H e [Hi,H2] 
with Hi > 0. The subject of this paper is the analysis 
of the metamagnetic response in correlated electron sys- 
tems in terms of the Fermi liquid description JT]). For 
this we calculate the effective mass and the term due 
to quasiparticle interactions from a microscopic model. 
This allows us to understand what drives the magnetic 
response. This can be relevant for the interpretation of 
experiments for itinerant metamagnets where the mag- 
netic response is measured simultaneously with the field 
dependence of the specific heat. 

In a naive single electron picture itinerant metamag- 
netism is not intuitive as with increasing polarization 
the magnetic response usually decreases. For instance, 
in weakly interacting systems, such as a Hubbard model 
with small U, with a featureless concave density of states 
metamagnetic behavior does not occur. RPA based cal- 
culations yield a decreasing susceptibility with increasing 
field as spin fiuctuations are suppressed. On the other 



hand, a convex density of states, i.e. with positive cur- 
vature at the Fermi energy, such as in the Wohlfahrt 
and Rhodes^ theory can lead to metamagnetic behav- 
ior. This is exploited in a number of works, where the 
Hubbard model with such convex density of states is 
analyzed-:-. Metamagnetic behavior is shown to also oc- 
cur in situations where the Fermi energy lies close to a van 
Hove singularity^!^, or where a Pomeranchuk Fermi sur- 
face deformation instability occurs^. It has been shown 
by calculations based on the Gutzwiller approximation 
by Vollhardlii and Spalek and coworker a^i^i^° that for a 
generic concave density of states metamagnetic behav- 
ior is also found in the intermediate coupling regime of 
the Hubbard model. The metamagnetic scenario is then 
that of correlated electrons, with a (Mott) localization 
tendency due to the interaction. 

Our calculations are based on the half filled sin- 
gle band Hubbard model which has been used fre- 
quently to describe itinerant metamagnetism for cor- 
related electrons^i^i^i^i^'iiii^ due to its relative formal 
simplicity. We employ the dynamical mean field the- 
ory (DMFT) ^^:^^ combined with the numerical renor- 
malization group (NRG}^'i^ to solve the effective impu- 
rity problem. We focus on the case of zero temperature, 
where sharp features are most clearly visible. We follow 
these earlier approaches here and restrict ourselves to the 
response of the paramagnetic solutions of the Hubbard 
model, which is possible for mean field-like approaches. 

The half filled Hubbard model in a magnetic field has 
already been investigated by detailed DMFT studies by 
Laloux et al.— and Bauer and Hewsooi^. Low tem- 
perature magnetization curves and field induced metal 
insulator transitions have been investigated by Laloux 
et al. Metamagnetic response based on correlated elec- 
tron physics, seen in the Gutzwiller approach, was con- 
firmed in such calculations. Our analysis extends previ- 
ous workii as we investigate the T = magnetic response 
with a Fermi Hquid interpretation based on the field 
dependent renormalized parameter approachi^^ii^ii^. 



This, together with results for the spectral functions, al- 
lows us to identify what gives rise to the magnetic re- 
sponse in the system. 

The paper is organized as follows. In a brief section II 
we give details about the model and method. The Fermi 
liquid interpretation and the relation between Fermi liq- 
uid parameters and the field dependent renormalized pa- 
rameters are described in section III. Section IV reports 
the results for magnetization, susceptibilities and the in- 
terpretation in terms of effective mass and quasiparticle 
interactions. We conclude by putting our results in con- 
text with itinerant metamagnetism studied experimen- 
tally. 



II. MODEL AND METHOD 

The basis for our calculation forms the Hubbard 
Hamiltonian in a magnetic field, which in the grand 
canonical formulation reads 



^M = Xl(*y4'^'^j> + ^■*^-) - Xl^^"*" 



u 



^»^vT"j,i- 



i,J,o' 



(2) 
c]^ creates an electron at site i with spin cr, and rii^a — 

cl^Ci^cr- tij = —t for nearest neighbors is the hopping 
amplitude and U is the on-site interaction; ^o- = M + ^^^ 
where fi is the chemical potential of the interacting sys- 
tem, and the Zeeman splitting term with external mag- 
netic field H is given hy h — g^^H/2 with the Bohr 
magneton /iB. In the DMFT approach the proper self- 
energy is a function of lo only^fiiSi. In this case the local 
lattice Green's function G'^°'^{io) can be expressed in the 
form, 



GTiu^) = de 



Po(e) 



Mc 



-S^H-e' 



(3) 



where po{s) is the density of states for the non- interacting 
model ([/ = 0). It is possible to convert this lattice prob- 
lem into an effective impurity one^^, introduce the dy- 
namical Weiss field QqI{lj). The DMFT self-consistency 
condition reads 



c?o:^H = GirM-^ + s.M. 



(4) 



The Green's function G^°^{iu) can be identified with the 
Green's function Ga{to) of an effective Anderson model, 
and Qq]^{uj) expressed as 



Q^» 



+ ^ia - K„{uj). 



(5) 



The function Ka{ijj) plays the role of a dynamical mean 
field describing the effective medium surrounding the im- 
purity. K„{ijj) and Sct(w) have to be calculated self- 
consistently using equations |[3l)-|[l]). Our calculations 
are based on the numerical NRG^i^^ to solve the effec- 
tive impurity problem. As in earlier worki^ we calculate 



spectral functions from a complete basis se t^^i^^ and use 
higher Green's functions to obtain the self-energy2i. For 
numerical calculations within the DMFT-NRG approach 
for po(£) we take the semi-elliptical form for the non- 
interacting density of states Po°™(^) = SV^*^ — e^/ttZ?^, 
where W — 2D is the band width with D — 2t for the 
Hubbard model, t — I sets the energy scale in the fol- 
lowing. 



III. FIELD DEPENDENT RENORMALIZED 
PARAMETERS AND FERMI LIQUID THEORY 

The response of a metallic system of correlated elec- 
trons can often be described in terms of Fermi liquid 
theory. The ratio of the spin susceptibility of the inter- 
acting system Xs and that of the non-interacting system 
Xs is given in equation ^. Thus, when strongly inter- 
acting fermions have a large paramagnetic susceptibility, 
it can be interpreted as due to quasiparticles with large 
effective masses. It is, however, also possible that the sus- 
ceptibility is additionally enhanced due to the quasipar- 
ticle interaction term 1/[1-|-Fq], which is for instance the 
case in liquid "^He, where m* /mo ~ 5 but Xs/Xs — 20?^ 
This is usually described by the dimensionless Sommer- 
feld or Wilson ratio R of the magnetic susceptibility and 
the linear specific heat coefficient 7. We will use it in the 
form R = (Xs/x°)/(7/7o), where 7/70 = m*/mo. 

Here we are interested in analyzing the behavior in 
finite field, and it is possible to calculate corrections of 
higher order in H to equation JTJjf^ We will, however, fol- 
low a different approach here, and assume that expression 
(U) remains valid for finite field with field dependent ef- 
fective mass m*{H) and Landau parameter Fq{H). This 
is in the spirit of the field dependent quasiparticle param- 
eters introduced in earlier worki^^ii^. Notice that for 
the case considered the field dependence of Xs, which is 
given by the non-interacting density of states, varies very 
little in the relevant field range. In this picture with field 
dependent parameters, metamagnetism can occur when 
the effective mass increases with the magnetic field. Gen- 
erally, however, also the field dependence of the quasipar- 
ticle interaction plays a role. One hypothesis, tested in 
this paper, is that itinerant metamagnetic behavior is al- 
ways accompanied by a field induced localization and a 
sharp increase of the effective mass near the metamag- 
netic transition. 

In order to calculate the microscopic Fermi liquid pa- 
rameters, we expand S(t(w) in powers of to for small uj, 
and retain terms to first order in w only. This is used to 
define renormalized parameters'^ 

fio,a = ZaW - ^am, and z, = 1/ [1 - E^ (0)] . (6) 

and from |(3]) a normalized quasiparticle propagator. 



GZi^) 



de 



' + flo,a - £ 



(7) 



Note that this w-expansion can also be carried out in 
finite magnetic field. Then the renormalized parameters 
become field dependent, Za — Za{h) and /io,(T = P'0,a{h). 
The density of states po,o-(e) derived from ([7]), po.ais) = 
—lTaGo^cr{£+'iS)/TT = po[{e+p,o^a)/ Za]/ Za, is referred to as 
the free quasiparticle density of states. Za- is interpreted 
as the weight of the quasiparticle resonance and /io,^ gives 
the position of the quasiparticle band. All energies are 
measured from the chemical potential fi. 

To obtain the renormalized parameters Za- and /io.o- , we 
use two different methods based on the NRG approach. 
The first method is a direct one where we use the self- 
energy S(j(w) determined by NRG and the chemical po- 
tential fia-, and then substitute into equation ([6]) for z^ 
and /io.cr- The second method is indirect, and it is based 
on the quasiparticle interpretation of the NRG low energy 
fixed point of the effective impurity.— This approach has 
been used earlier for the Hubbard modeU^i^ and for the 
Anderson impurity model in a magnetic fiel d^^'^^ . As 
shown before the results of both methods usually agree 
within a few percent, and we use an average value of both 
methods for the numerical results presented later. It is 
important to calculate these parameters accurately, since 
for the following results also their derivatives are needed. 

We can calculate static expectation values and re- 
sponse functions in terms of the renormalized parame- 
ters. The quasiparticle occupation number n° is given 
by integrating the quasiparticle density of states up to 
the Fermi level, 



n' 



= / de po,o-(e) == de po,cr{£)d{ncr - ^a 



s). (8) 



Luttinger's theorem^^ holds for each spin component for 
the Hubbard model in magnetic fieldri^, hence we have 
n^ = no-, where Ua- is the value of the occupation number 
in the interacting system at T = 0. 

To calculate the magnetic response we focus for the 
rest of this paper on the case with particle-hole symmetry 
where /i = U/2, and we can write Scr(0, h) ^ U/2—ar/{h). 
We can calculate j]{h) directly from the self-energy, e.g. 
ri{h) = (S| — S|)/2, or from the renormalized parameters 
ri{h) — p,Q{h)/z{h) — h. At half filling we have z^ ~ zi = z 
and /io.T = ^fk),i = Ao- We define the function 

g{h) := h + r]{h) = jiii{h)/z{h) = p,o{h)m*{h)/mo, (9) 

as m* /mo — z^^ in DMFT. In terms of the quasiparti- 
cles it is the product of the effective mass enhancement 
TO* /mo and the shift of the quasiparticle band /iQ. With 
the applicability of Luttinger's theorem the magnetiza- 
tion is then given by 

oo 

ni{h) = \{n^- n^) = I As Po{e)e\g(K) -e\-]^. (10) 



For a local self-energy this is an exact expression for the 
magnetization, which only depends on the field depen- 
dent renormalized parameters via g(}i). For certain bare 



densities of state, for instance, for the semi-elliptical den- 
sity of states p^Q^ie), it can be evaluated analytically, 

m{h) - \g{h)p^r{9{h)) + - arcsin(g(/i)). (11) 

Differentiating lfTO|) with respect to h yields the local 
static spin susceptibility 



Xs = -^ = g'{h)po{g{h)) 



(12) 



where here and in the following primes indicate deriva- 
tives with respect to h. A similar expression had already 
been derived by Luttinger— . The metamagnetic condi- 
tion x'sW > is then 

g"ih)poigih)) + p'Mh))9'{hf > 0. (13) 

The occurrence of metamagnetic behavior can be ana- 
lyzed depending on the functional form oi g{h) and pa^e). 
For a simple analysis let us assume h > and the power 
law form for g{h) = ch", c > 0. The first term in lfT3|) 
is then positive if a > 1. For a convex density of states, 
Po{s) > 0, the second term is also positive and meta- 
magnetic behavior occurs as mentioned earlier. For a 
concave density of states, Po(e) < 0, the two terms in 
(fT3l) compete. If we also assume the power law form for 
the density of states, po{£) ~ r^ — de^ , (e.g. for p5°™ 
one has tq = 2/itD d — ^0/2 and 7 = 2) condition lfT3|) 
becomes 



ro 



1 



c'>'da(l+7)- 1 



>/i"^ 



(14) 



Since the right hand side is positive, we can infer that for 
a > 1 and 7 > (1 — a) /a metamagnetic behavior occurs. 
The actual field dependence of g{h) can be calculated 
from the renormalized parameters and it depends on the 
interaction strength. As we will see for the half filled 
Hubbard model and intermediate U , g{h) grows faster 
than linear with h, i.e. a > 1. 

In the limit of zero field the ratio of the susceptibility of 
the interacting and non-interacting system has a simpli- 
fied expression in terms of the renormalized parameters. 



XS ,,„^ "i*(0) , 

— =5(0) = —z — Aio(O): 



X 



Too 



(15) 



for /io(0) = 0. Comparing with the Fermi liquid expres- 
sion (U) we can identify 1/(1 -I- Fq) = Pq. This quantity 
corresponds to the Wilson ratio R. In the general case, 
the field dependent enhancement due to the quasiparticle 
interactions reads 



R{h) 



1 



l + FSih) 






m^ 



Poih) 



So far the considerations have been independent of 
our DMFT-NRG approach. In the following section we 



will compare results for the magnetic susceptibility ob- 
tained from the static expectation values of integrating 
the Green's functions, with the results based on the field 
dependent parameters. We determine them as described 
above. Alternatively they can be calculated by other 
methods, such as the Gutzwiller (GW) approach, and we 
will make comparison as appropriate. Results are ob- 
tained as in Ref. [tI, where the critical interaction for the 
metal insulator transition is U^^ = 16W/3tt « 6.79 for 
p^™(e) with W = 4. 



IV. RESULTS 

A. Magnetization and metamagnetic transition 

For a first overview we present results for the magne- 
tization m{h) as a function of field h in Fig. [T]for vari- 
ous values of U. The magnetization m{h) was computed 
from the static NRG expectation value (EV) for the oc- 
cupation number as well as from integrating the spectral 
function to the Fermi level, both of which agree very well. 
The results for m{h) based on the field dependent renor- 
malized parameters (RP) and equation ifTTj) are also in 
good agreement, but not included in the figure. 




FIG. 1: (Color online) The local magnetization m{h) as a 
function of the magnetic field h for different values of U. We 
can see that a metamagnetic curvature sets in at (7 = 3. Inset: 
Hysteresis curve for U = 4 (triangle up increasing h, triangle 
down decreasing h). 

The plot gives a clear picture of the field strength /ipoi 
necessary to polarize the metal completely to ?ti = 1/2. 
For weak coupling it can be related to the rigid band 
shift and a large field h ^ D is needed, but for larger 
interaction strength /ipoi is reduced substantially. For 
U > 3 a metamagnetic curvature in the magnetiza- 
tion can be observed, and we see that in the Hubbard 
model at zero temperature the metamagnetic transition 
field2S ft-m coincides with /ipoi, which is not necessarily 



the case for T > 0. Laloux et edM- have compared re- 
sults from low temperature DMFT calculations with the 
Gutzwiller approximation and it was found that the oc- 
currence of metamagnetic behavior is overestimated by 
the Gutzwiller approximation (see also Fig. [3]). 

Earlier workii showed that the transition is a discon- 
tinuous first order one at low temperature. Our results 
show jumps in the magnetization curve at the transition 
field ft-m, e.g. for U — 3 and U — Am Fig. [H however, we 
can not exclude a very steep continuous increase which 
can not be resolved numerically. We have also found hys- 
teresis, shown for C/ = 4 as an inset in Fig. [T] (triangle up 
increasing h, triangle down decreasing h). This suggests 
that the transition is also of first order for zero tempera- 
ture. For larger interaction U > 4.5 there exists a small 
field range near /im, where we have not found unique, well 
converged DMFT solutions, so no definite statement can 
be made. 

The half filled repulsive Hubbard model in magnetic 
field can be mapped to the attractive one^S, in which the 
chemical potential is related to the field in the original 
model, fi = U/2 + h. The attractive model has been stud- 
ied by the DMFT in situations, where superconducting 
order was not allowed fo r^^'^^ . A first order transition 
from a metaUic to a pairing state for fixed density was 
found at a critical interaction. The occurrence of the 
transition can be related to the metamagnetic transition 
here. A nearly polarized system corresponds to a low 
density limit, and to estimate when the transition sets 
in, one can analyze the two-body problem in the attrac- 
tive model and calculate the critical Uc for bound state 
formation. For a three dimensional cubic lattice the re- 
sult is Uc « 0.659W^. With the given bandwidth this 
corresponds to a value of Uc ~ 2.64, which is a reasonable 
estimate for the interaction strengths, where the metam- 
agnetic behavior is found here. 



B. Magnetic susceptibilities and quasiparticle 
properties 

From the initial slope of the magnetization curves in 
Fig. [1] we observe an increase of the magnetic suscepti- 
bility with the interaction strength U. This increase can 
also be seen in the following Fig. [2] where we show the ra- 
tio of zero field susceptibility to the non-interacting value 
X° as function of U deduced from differentiating the EV 
for m{h) in the limit h ^ 0. 

For comparison we have also included the susceptibil- 
ity calculated from equation lfT5|) with the renormalized 
parameters (RP) and their derivatives, as well as the re- 
sults obtained from the Gutzwiller (GW) approximation. 
EV and RP results agree very well, confirming the appli- 
cability of Fermi Hquid results in this metallic regime. 
The GW results follow a similar trend but overestimate 
the value for the susceptibility, which becomes more pro- 
nounced for larger U. 

The inset plot shows the [/-dependence of the effective 
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FIG. 2: (Color online) The [/-dependence of the magnetic sus- 
ceptibility Xs- We compare results deduced from the EV of 
m{h) with ones obtained from the RP and from the Gutzwiller 
(GW) approximation. The inset shows the effective mass 
m*{U)/mo and the Wilson ratio R{U) as a function of U. 



mass and the Wilson ratio. In terms of Fermi liquid the- 
ory and the expression (U) the increase of Xs with U can 
be understood by the behavior of the effective mass and 
the progressive localization tendency, which brings out 
more the spin degrees of freedom of the electrons. We 
can see, however, that the effective mass ratio is larger 
than that of the magnetic susceptibility. This difference 
can be attributed to the factor R = JIq = [1 + Fq]~^, 
which is due to the quasiparticle interaction. This factor 
is larger than one for smaller values of U, but decreases 
to values below one for stronger interaction. This indi- 
cates a sign change of the parameter Fg from negative 
to positive. The comparison of the corresponding quan- 
tities calculated in the GW approximation shows a qual- 
itatively similar behavior for both m* /rriQ and R, when 
U is small. For larger values of U in Fig. [21 however, 
the effective mass enhancement in the GW approach, 
m* /mo = 1 — (U/U^^)'^, is much smaller and R in- 
creases with U in contrast to the DMFT result. 

We return the finite field response and focus on the 
metamagnetic behavior which is found for intermediate 
values of U. Results for the ratio of the magnetic suscep- 
tibility in finite and zero field deduced from differentiat- 
ing the magnetization (EV) are compared to the ones ob- 
tained from the quasiparticle parameters (RP) and equa- 
tion lfT2|) . For completeness, we have also included results 
from the GW approximation. This is shown in Fig. [3]for 
U — 3 m the upper panel and U — 4.5 in the lower panel. 
We can see that also in finite field the results for the sus- 
ceptibility calculated from the EV for m(/i) and the field 
dependent RP agree fairly well with a deviation of less 
than 5%. For the case U = 3 (upper panel) the results for 
x{h) based on the field dependent RP are always smaller. 
In both cases we find first a period where the susceptibil- 
ity is nearly constant, but then starts to increase rapidly 
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FIG. 3: (Color online) The /i-dependence of the ratio of the 
finite and zero field magnetic susceptibility Xs for U = 3 (up- 
per panel) and U = 4.5 (lower panel). We compare results 
deduced from the EV for m{h) with ones obtained from the 
RP and the ones from the GW approach. The inset shows 
the ratio of finite and zero field effective mass m*{h)/mo{0) 
and the Wilson ratio R{h)/R{0) as a function of h. 



as h approaches h^i. For U = 3 the values obtained from 
the RP initially decrease slightly with the field, which is 
however incorrect, and comes about through numerical 
inaccuracies when determining the parameters and the 
numerical differentiation. As /ij„ = /ipoi the magnetic 
susceptibility is zero for h > h^. At finite temperature 
a susceptibility maximum is expected. The results for 
Xs from the GW approximation show generally a similar 
trend, but as mentioned earlier the metamagnetic behav- 
ior sets in at lower field strengths. 

A difference in the behavior between the two cases is 
visible in the two insets where the ratios of field depen- 
dent effective masses to their zero field values and the 
field dependent Wilson ratios R{h)/R{0) are plotted. For 
the C/ = 3 case the effective mass decreases with the field 
which is typical behavior in the weak coupling regime. It 
can be understood by RPA approximations where spin 



fluctuations, which give an effective mass enhancement, 
are suppressed in finite field. The metamagnetic increase 
of the susceptibiUty, however, can not be explained by 
this. In terms of Fermi liquid theory it is related to the 
magnetic field dependence of the quasiparticle interac- 
tion rather than the localization tendency encoded in the 
effective mass. R{h)/R{0) indeed is increasing sharply 
close to /im- In equation (fT6|) we have two competing 
terms for this enhancement factor, m*' /m* < 0, but one 
finds JI'q > |/io™*'/"^*l which leads to the observed en- 
hancement. The drive for the metamagnetic behavior is 
therefore due to the shift of the quasiparticle band from 
the Fermi level with increasing field. This contrasts to 
the weak coupling situation, such as C/ = 2, where R{h) 
decreases with the field strength and no metamagnetic 
response is observed. 

The effective mass in the case oi U = 4.5 (lower panel 
in Fig. [3]) shows different behavior. We can see a sharp 
increase with the field. However, the magnitude the ra- 
tio m* /ruQ increases is less than that of the susceptibil- 
ity. The difference again can be related to the Fermi 
liquid factor R — 1/[1 -I- Fq], which is larger than one 
and increasing with h as can be seen in the inset of the 
lower panel in Fig. [3l In this case the second term in 
equation (flGl) is positive and the first term negative, but 
lAol < Imo"i*'/"^*I- The results from GW approach for 
the effective mass and R are in line with the DMFT cal- 
culations for the case U — 3, however, for U — 4.5, the 
GW result for m*' /m* only increases very little with the 
field, whereas R{h) increases sharply to yield the meta- 
magnetic response. 

For larger interactions than the ones discussed here 
(5 < C/ < C/c), one can encounter difficulties to reach 
convergency in the DMFT calculations with finite field 
as discussed in earlier worki^. The results indicate, how- 
ever, that there is a strong field dependent enhancement 
of the effective mass which is the main drive for the meta- 
magnetic response. The ratio R{h)/R{0) varies little with 
h or even decrease for larger fields. Such a behavior is 
also found within the GW approach for larger U near the 
metal insulator transition. 



C. Spectral functions 

The behavior of the quasiparticle band can be seen 
directly in the local spectral function. For the cases with 
smaller coupling the field dependent response shows a 
continuous shift of spectral weight to lower energies for 
the majority spin (see Fig. [4] for U = 2). 
Note that the minority spin density of states Pi{uj) is 
given by p^{—uj) at half filling. To illustrate the behavior 
of the quasiparticle peak for the stronger interacting case 
with U — 4.5 in more detail, we plot the local spectral 
function for the majority spin Pt('-^) i^ Fig. [H 
In the upper panel we can see how the lower Hubbard 
peak in the spectral density acquires weight when the 
field and thence magnetization is increased whilst the 




FIG. 4: (Color online) The majority spin density of states for 
(7 = 2 and various field strengths in comparison. 
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FIG. 5: (Color online) The majority spin density of states 
for U = 4.5 and various field strengths in comparison: upper 
panel full frequency range, lower panel low frequency behav- 



upper Hubbard peak loses spectral weight. The behavior 
at low energy is seen more clearly in the lower panel. At 
first sight the overall picture is reminiscent of the particle 
hole symmetric Anderson impurity model in the Kondo 
regime in magnetic fieldi^ as far as the high energy be- 
havior is concerned. The quasiparticle resonance in the 
locally correlated system broadens and departs from the 
Fermi level. This behavior occurs in an analogous fash- 
ion in the weak coupling regime of the Hubbard model 
with /io(/i) > 0. In the strongly correlated case, however, 
we find a significant narrowing of the quasiparticle peak 
in the field, which is accompanied by the field induced 
metal insulator transition and metamagnetic behavior. 
The quasiparticle resonance first departs from the Fermi 
energy, but for larger fields is driven back to it. These 
features are visible in the field dependence of the renor- 
malized parameter /io with JIq < as discussed above. 



V. RELATION TO EXPERIMENTS AND 
CONCLUSIONS 

It is of interest to see, whether the described be- 
havior bears any resemblance with what is observed 
experimentally in strongly correlated itinerant electron 
system. Metamagnetic behavior is observed, for in- 
stance, in the heavy fermion compounds CeRu2Si2^'^, 
UPts^ or Sr3Ru207^'^'22i38 and the Co-based metallic 
compounds such as Y(Coi_a;Ala;)2,— >^ sometimes called 
nearly ferromagnetic metals. The microscopic origin for 
the occurrence of the effect in these compounds can be 
manifold, and is sometimes still controversial. In many 
cases antiferromagnetic exchange is thought be impor- 
tant and the system's closeness to a magnetic instability. 

For generic features, we attempt to compare our micro- 
scopic Fermi liquid description with experimental stud- 
ies of itinerant metamagnetic behavior in heavy fermion 
compounds. It is important, however, to be aware that 
our results based on the paramagnetic solutions of the 
half filled single band Hubbard model are not appro- 
priate to make quantitative predictions for those com- 
plex systems. Organic conductors are thought behave 
like simple Mott-Hubbard systems and have been shown 
to display a magnetic field induced localization transi- 
tion with hysteresis by resistance measurements.— The 
author is, however, not aware of any published field de- 
pendent magnetization or specific heat data to compare 
to. 

In materials such as CeRu2Si2, UPta or Sr3Ru207 
the magnetic field dependence of the linear specific 
heat coefficient 7 was measured near the metamagnetic 
transition^i^i^^i^. It is worth noting that, as can be 
shown from a thermodynamic identity, the field depen- 
dence of 7 can also be extracted from T^-coefHcient of the 
magnetization^. In the experiments 7 increases with the 
magnetic field and possesses a maximum at the metam- 
agnetic transition h = h^. This is comparable with the 
Fermi liquid results for stronger coupling, e.g. the case 



U = 4.5 (Fig. [3]lower panel), where the effective mass in- 
creases with the magnetic field. In the case of CeRu2Si?2i 
one can see that the susceptibility increase with the mag- 
netic field is up to about 8.5 times the zero field value, 
whereas in the same regime the specific heat coefficient 
only shows an enhancement of 1.6. In our Fermi liq- 
uid interpretation this signals that the quasiparticle in- 
teraction plays an important role in the susceptibility 
enhancement and the metamagnetic behavior. The rel- 
evance of this has been emphasized in the recent exper- 
imental work on Yb3Pt4.~ A more careful quantitative 
comparison would be possible based on the periodic An- 
derson model, for instance. The presented approach can 
be extended to this situation, but also other techniques 
are available^2iiid^ii£. 

To summarize, we have analyzed the metamagnetic 
response of the half filled Hubbard model in terms of 
renormalized quasiparticle parameters and Fermi liquid 
theory. The renormalized parameters can be calculated 
accurately with methods based on the NRG, and they 
have a clear physical meaning. It is shown that the field 
dependent metamagnetic behavior can have part of its 
origin in field induced effective mass enhancements, but 
is not fully explained by this. This is most clearly pointed 
out in Fig. [Sj where metamagnetic behavior for smaller 
U is accompanied by an effective mass reduction in the 
field, whereas for larger interaction the opposite is the 
case. The comparison with results obtained from the 
Gutzwiller approximation gives similar trends, but shows 
quantitative deviations. The hypothesis that the meta- 
magnetic behavior in itinerant systems is always driven 
by field induced mass enhancement is therefore found to 
be not valid. In the intermediate coupling regime it is 
also shown that the effective mass enhancement alone is 
not sufficient to explain the metamagnetic enhancement 
and based on Fermi liquid theory arguments the quasi- 
particle interaction has to account for the difference. As 
a generic feature there the corresponding term described 
by the Wilson ratio R increases near the metamagnetic 
transition. The opposite happens in the weak (no meta- 
magnetic response) and strong coupling situation. The 
observation that only a part of the susceptibility enhance- 
ment is based on the effective mass is found to be qual- 
itatively in agreement with experimental observations in 
heavy fermion systems. 
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